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Abstract: We investigate the interactions of large composite dark matter (DM) states 
with the Standard Model (SM) sector. Elastic scattering with SM nuclei can be coherently 
enhanced by factors as large as where A is the number of constituents in the composite 
state (there exist models in which DM states of very large A > 10® may be realised). 
This enhancement, for a given direct detection event rate, weakens the expected signals 
at colliders by up to 1/A. Moreover, the spatially extended nature of the DM states 
leads to an additional, characteristic, form factor modifying the momentum dependence 
of scattering processes, altering the recoil energy spectra in direct detection experiments. 
In particular, energy recoil spectra with peaks and troughs are possible, and such features 
could be conhrmed with only 0(50) events, independently of the assumed halo velocity 
distribution. Large composite states also generically give rise to low-energy collective 
excitations potentially relevant to direct detection and indirect detection phenomenology. 
We compute the form factor for a generic class of such excitations—quantised surface 
modes—finding that they can lead to coherently-enhanced, but generally sub-dominant, 
inelastic scattering in direct detection experiments. Finally, we study the modifications 
to capture rates in astrophysical objects that follow from the elastic form factor, as well 
as the effects of inelastic interactions between DM states once captured. We argue that 
inelastic interactions may lead to the DM collapsing to a dense configuration at the centre 
of the object. 
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1 Introduction 

Most models of dark matter (DM) assume that it can, for practical purposes, be treated as 
a collection of point-like particles. However, this is not necessarily the case, and a variety of 
models in which DM is a composite state have been proposed (for example, WIMPonium 
[1-3], and dark atoms [4]). Given the centrality of DM to our present thoughts about 
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beyond-the-Standard-Model physics, and our present lack of knowledge concerning many 
of its fundamental characteristics, and thus how DM may reveal itself in experiments or 
observations, it is important to consider possible variations away from the standard picture. 
This is especially the case if there are qualitatively new features compared to traditional 
DM models that may affect direct and/or indirect detection phenomenology. 

In this paper, we consider the consequences for present-day scattering processes of a 
simple kind of compositeness, in which the DM states are composed of a large number. A, 
of constituents, forming an extended semi-uniform object, analogous to Standard Model 
(SM) nuclei formed out of constituent nucleons.^ Models that realise such a scenario 
include Q-balls (non-topological solitons carrying a conserved charge) [6, 7], and Nuclear 
Dark Matter models, in which DM is made up of bound states of strongly-interacting 
constituents with short-range interactions [8-10]. The early-universe cosmology of these 
models has interesting features, as investigated in [7-10]; here, we simply assume that a 
late-time population of such states exists. 

This large and extended compositeness affects elastic scattering and the associated 
direct detection phenomenology, as we discuss in Section 2. Large composite states con¬ 
taining A constituents can have low-momentum-transfer elastic scattering cross sections 
coherently enhanced by a factor as large as A^ [8, 10, 11]. This case is realised if the size of 
the state is not too large compared to the inverse momentum exchange relevant in direct 
detection scatterings. In this situation, and assuming that the mass of the composite state 
is oc A, so the DM number density is oc 1 /A, the event rate at direct detection experiments 
will effectively be enhanced oc A for a given interaction strength between SM and DM 
constituents. Previous studies of the build-up of composite DM states in the early uni¬ 
verse show that large values of A (> 10®) could plausibly be realised [10], so the effective 
enhancement can be significant. 

In addition, if a high enough proportion of the DM states have radii larger than SM 
nuclei, but not so large as to significantly suppress coherent scattering, the spatial extension 
of the DM states leads to a dark form factor modifying the momentum dependence of 
DM-SM scattering (as previously considered in [11]). In the simplest case of a scalar 
interaction depending only on the density, this form factor has a characteristic series of 
peaks and troughs, analogous to SM nuclear form factors. In direct detection experiments 
with sufficiently good energy resolution, these could lead to the striking signature of rises 
in the energy recoil spectrum. We find that significant features of this kind could be 
distinguished from point-like elastic scattering after the observation of only 0(50) events, 
independently of assumptions about the DM halo velocity distribution. A distribution 
over DM sizes may average out these peaks into a smoothly-falling effective form factor, 
in many cases resembling that from e.g. the exchange of an intermediate-mass mediator, 
but multiple direct detection experiments using different SM nuclear targets would still be 
able to separate out the momentum dependence from halo velocity distribution effects, as 
with other models of dark form factors [11-15]. 

^ The generic scattering phenomenology of small-number composite DM states has been considered 
in [5]. 
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As we discuss in Section 3, large composite states also generically give rise to long- 
wavelength and low-energy collective excitations. There will also be higher-energy exci¬ 
tations, and fission processes, but making the binding energy scale low enough for these 
to be relevant in direct detection experiments generally complicates the cosmology of such 
models [10]. We focus on the most model-independent possibility of collective density ex¬ 
citations. We find that they can be of sufficiently low energy to be excited in collisions 
with SM nuclei, leading to coherently-enhanced inelastic scattering. Although this is sub¬ 
dominant to elastic scattering in the scenarios we consider, a non-negligible fraction of 
events may be inelastic, and we calculate the associated form factor in leading approxima¬ 
tion.^ 

As is well known, couplings to SM nuclei may also lead to the capture of DM by 
astrophysical objects. As discussed in Section 4 we find that the composite nature of the 
DM changes the overall capture rate. In addition, perhaps the most interesting qualitative 
effect is that the large local DM density inside the star, and the natural possibility of 
inelastic DM self-interactions, could lead to a DM distribution that is significantly modified 
compared to the naive expectation. In particular, for either dissipative collisions or fusions, 
there is the possibility that this process runs away to a state in which most of the captured 
DM lies in a single very dense conhguration.^ There may be model-dependent consequences 
of the energy released by inelastic DM self-interactions, which in the case of fusions could 
release large amounts of energy into small volumes over short timescales. 

2 Modifications to elastic scattering and direct detection 
2.1 Dark sector form factors 

For elastic scattering between a point-like state and a spatially extended state, interacting 
via short-range interactions sourced by some density p on the extended side, the dependence 
on the spatial properties is summarised by the form factor 

F(q) = I dre*^-p(r), (2.1) 

i.e. by the Fourier transform of the density. For scattering of DM states off SM nuclei, 
effective field theory arguments show that the scattering operator on the SM side should 
take on one of a restricted number of forms (see e.g. [16] for a recent comprehensive anal¬ 
ysis). The dark sector density must have a complementary tensorial structure, but could 
in general be determined by any properties of the state, e.g. density, spin, etc. 

^Taking the possibility of inelastic DM seZ/-interactions further, we may also expect, in the context of 
galactic halos, exothermic DM-DM interactions, e.g, from fusions. This can lead to an increase in the 
average kinetic energy of DM in the galactic halo, which might in turn result in a clearing out of the central 
high-density region. This possibility was outlined in [10], and requires a dedicated study of halo dynamics. 
In this work, we focus on the interactions of composite DM states with the SM, and leave the topic of 
modified halo dynamics to a future paper. 

®Since the large composite states with saturated densities we consider in this paper require a short-range 
repulsive interaction between constituents, forming a black hole at the centre of the star would generally 
require accumulating a very large DM mass. 
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In the case of two spatially extended states scattering off each other, the effective 
density is the convolution of the separate spatial profiles, so the matrix element is found by 
multiplying the form factors together. Thus, the overall form factor for DM-SM scattering 
will be FAr(q)Fx(q), where Fn and Fx are, respectively, the SM nuclear form factor and the 
dark sector form factor (as we will discuss below, in many experimental circumstances the 
SM nuclear form factor is of limited importance). The specihc density determining the form 
factor will depend on the nature of the DM-SM interaction. For simplicity, in this paper 
we will restrict ourselves to considering the case of a scalar form factor depending only on 
the number density of the state, e.g. arising from the exchange of a heavy scalar mediator 
coupling uniformly to all of the constituents. Taking, as discussed in the Introduction, the 
the composite DM states to consist of an approximately uniform density of constituent 
matter, a hrst approximation for the density is a spherical top hat function, leading to a 
spherical Bessel function form factor, 


^ SAMqR) 
qR 


3A{sm{qR) — qRcos{qR)) 


( 2 . 2 ) 


where R is the radius of the top hat density, and A is the total volume integral of the 
source across the distribution, so F{0) = A. 

SM nuclei provide an example of this kind of roughly-constant-density state, and illus¬ 
trate the kind of deviations from the top-hat form factor that might occur. Figure 1 shows 
an example of the form factor corresponding to the nuclear charge density for the particular 
isotope ^°Ge [17], as inferred from electron scattering data. Generally, and in this specific 
case, the hrst few peaks and troughs of SM nuclear form factors are well-approximated by 
the ‘Helm’ functional form 


^ An{qR) 

qR 
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which is simply the top hat form factor modihed to have a hnite-width fall-off, over the 
‘skin depth ’s ~ 0.9fm (comparable to the scale of the individual nucleons). 

Importantly, in the case of interest to us where both a dark-sector and SM form-factor 
are present, if the dark states have larger radii than SM nuclei, then the hrst zeros, or more 
generally troughs, of Fx{q) will occur at smaller q than for Fx{q)- This means that, while 
the SM nuclear form factor is of limited importance in most direct detection experiments, 
the dark form factor may have interesting structure in precisely the momentum, and thus 
recoil energy, ranges being probed. 


2.2 Direct detection recap 


In this section, we briehy review the standard direct detection formalism leading to the 
differential event rate (equation 2.7), mainly to establish notation for what follows. Readers 
familiar with this material may safely proceed to Section 2.3. 

The differential scattering rate (event rate per unit target mass) at which incident DM 
particles scatter off an initially-stationary target SM nucleus, giving it recoil energy Er, is 
given by 
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da XN 
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Figure 1. Form factor for nuclear charge density (black) of ^°Ge nucleus, as derived from electron 
scattering data [17]. Blue (dashed) curve is standard Helm parameterisation, red (dot-dashed) 
curve is constant-density approximation. 


where v = jvj, rriM is the mass of the SM nucleus, nxf{^)d^'v is the differential number 
density of incident DM particles, Umin is the minimum velocity required to obtain recoil 
energy Eji, and (Txn is the DM-nucleus scattering cross section."^ For non-relativistic elastic 
scattering, Umin = where uxN = is the DM-nucleus reduced mass. 

The momentum transfer in an elastic collision is q = ii,xNV-\/2iX^-cosE^, where 9* 
is the scattering angle the CoM frame. Since Er = 5 ^^, we have (IEr = 
where dW is differential solid angle in the CoM frame. Referring back to the discussion of 
Section 2.1, the matrix element for scattering at angle 6* will depend on the momentum 
transfer q, and possibly also on the velocity v. Treating isotropic, velocity-independent 
scattering to start with, and writing -^crxNio) = we have 


dcrxN 

dEji 


mx 

‘^Exn'^'^ 


(yxN{q)- 


(2.5) 


So, under the assumption of velocity-independent scattering, we can factor the differential 
scattering rate as 


dR 

dEji 



/ \ _ / \ / \ 

-“2 - 0-XN[Q) = 






( 2 . 6 ) 


Altering notation somewhat from Section 2.1, we can write axNiq) = (^XNFN{q)‘^Fxiq)‘^, 
where axN is the zero-momentum-transfer cross section (including any coherence enhance¬ 
ment), so that |Fjv(0)| = |Fx(0)| = 1 . 

The scattering rate is usually expressed in terms of the DM scattering cross section 
with nucleons, instead of with full nuclei. For dimension -6 interactions between the DM 
constituents and SM quarks, the zero-momentum-transfer scattering cross section goes 

^Throughout we will consider parameter ranges where the DM-SM interactions are sufficiently weak that 
the Earth is optically thin to DM, so that nx/(v) is the DM halo distribution at the location of the Earth. 


- 5 - 























dR/dE,i dR/dEji 



Figure 2. Recoil energy spectrum (black, solid curves) for Bessel-function form factor, radius 50 fm, 
with DM velocity distribution as given by Standard Halo Model [18], assuming Left: Germanium 
detector target, with Gaussian energy response (cr(£’) = ^0.3^ + O.OO^Fl/keVkeV), and Right: 
Xenon detector target, with Gaussian energy response (cr(i?) = O.bkeVi/if/keV), assuming 5keV 
energy threshold for both. The DM state is taken to be composed of 3 x 10® constituents, each 
of mass 20GeV, with constituent-SM nucleon cross section of 2 x 10“^® pb. For comparison, blue 
(dashed) curves show energy recoil spectrum for a 20 GeV WIMP, red (dot-dashed) for a 1 TeV 
WIMP, both with axn = I0“® pb. Energy response functions are taken from projections for future 
experiments in [15]. 


2 

as where A is the suppression scale associated with the interaction, and C is the 

coherence enhancement factor. Here, this is given by the product of the integrated densities 
relevant to the microscopic interaction for the SM nucleus and the DM state, which we 
take to be simply the respective nucleon numbers. Thus axN = where uxn 

is the DM-nucleon reduced mass, and so^ 


dR 

dEji 




nx \Cn? 


fyxnFN{q)‘^Fx{q)‘^ ■ 


(2.7) 


2.3 Coherent enhancement of scattering rates 

With the assumptions made in the previous section, the zero-momentum transfer cross 

section between a DM state with A constituents and a SM nucleus with N constituents 
2 

is ao ~ Taking the mass of the DM state to be oc A, the DM number density 

is oc ^ oc assuming temporarily that all of the DM is of the same size. So, the 
overall scattering rate in direct detection experiments will, for mx ^ mx, he oc A for fixed 
SM-constituent interactions. This is true if the DM radius R is small enough that typical 
momentum transfers do not probe qR >1 — otherwise, the scattering rate is suppressed 
by the DM form factor, as discussed in Section 2.1. While, as discussed in Section 2.4, 

®This equation describes what a detector with perfect energy resolution would see. If a detector has 
some response function k such that the differential rate to detect events with true energy Eh, occurring 
at rate R, at measured energy E'h, is Rk[Eb., E'H)dE'ii, then the differential rate for measured events is 
RAEr) = f dE'HK{ER,E'H)^. 
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there might be some number distribution over DM states with different radii, if the mass 
distribution and scattering-rate distributions are confined to a small range in logarithmic 
size, then the oc A enhancement will hold approximately. 

Figure 2 shows an example of the direct detection recoil energy spectra resulting from a 
scenario of this kind, corresponding to DM particles with a Bessel-function form factor of ra¬ 
dius 50 fm. These are compared to the recoil spectra for standard momentum-independent 
(i.e. Fj{ = 1) WIMP scattering. Note that the constituent-SM nucleon cross section re¬ 
quired in the composite model is much smaller than that required in WIMP models giving 
approximately the same event rates. An additional point is that the much better energy 
resolution expected in solid-state experiments would enable them to resolve the peaks and 
troughs of a dark form factor corresponding to much larger radii than for liquid-phase (e.g. 
Xenon) experiments. 

The possible oc A enhancement means that, for a given direct detection event rate, the 
expected production of lighter single constituents in SM processes in colliders is reduced. 
For the example of DM coupling through the Higgs portal, the strongest collider constraint 
for DM states with mx < 'mh/2 comes from the Higgs invisible width. As described in [19], 
a bound of < 10% on the invisible branching ratio of the Higgs puts constraints on the 
DM-Higgs coupling a factor of a few better than current direct detection experiments for 
DM masses a small factor lower than m/i/2. Thus, only a modest relative suppression of 
collider rates vs direct detection rates is needed to render the direct detection bounds more 
constraining. Taking the example in Figure 2, with DM constituents of mass 20 GeV each 
having a nucleon scattering cross section = 2 x 10“^^ pb, we would require Higgs invisible 
width bounds of ~ 10“® to be competitive with current direct detection experiments. 


2.4 Recoil spectrum from size distribution 


The previous sections considered a dark form factor arising from DM particles having a 
single, common size. However, in many models of large composite states (e.g. [10], Q- 
balls), the cosmological process through which a population of these states arises generates 
a distribution over multiple sizes, which would lead to a smeared-out signature in energy 
recoil spectra. 

For DM consisting of a set of states Aj, the recoil spectrum for elastic scattering will 


be 


dR 

dEji 


FNiqf'^giiv min 
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Ui \Ci\‘^ 

2/4 


cri^nFx,iiqf, 


( 2 . 8 ) 


where the subscripted quantities are those of equation 2.7 for each species Aj. In the specific 
case of a spectrum of related states, this may simplify somewhat. Absent astrophysical self¬ 
interactions (see Section 4.2), all of the states may be expected to have the same velocity 

2 

distribution, i.e. gi = g. From cr„ ~ we expect the gi^n dependence to cancel, and 

additionally, if all of the states are heavy (m* 3> m^r), then — mx, so we can factor 
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out the g dependent term fully. Overall, 
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( 2 . 10 ) 


where nx = have replaced the sum over individual states with a single 

effective form factor, choosing C appropriately so that Fx{0) = 1. 

Figure 3 shows the effect of a distribution over DM states of the kind considered in [10]. 
The sum over different sizes corresponds to adding together differently scaled form factors, 
which smooths out the troughs of the Bessel-function form factor. In this case, we end up 
with a form factor which is similar to that obtained from the exchange of a light mediator 
particle, F{q) oc ; but arising entirely from contact interactions, without the need 

for a light state.® For a given DM velocity distribution, these kinds of dark form factor 
gives an energy recoil spectrum with a shape different from that of standard momentum- 
independent scattering, falling off more slowly at high energies than a low-mass WIMP, 
and more quickly at low energies than a higher-mass one. In the next section, we discuss 
how such form factors could be distinguished from momentum-independent scattering. 


2.5 Dependence on DM velocity distribntion 

From the definition in equation 2.6 of g'(uinin) = f ^ , it is clear that, by choosing 

/(v) appropriately, we can make g any non-increasing function of Umin (with the physical 
constraint that it must fall off very fast past Vg -|- Vesa where Vg is the velocity of Earth 
relative to the Galactic rest frame, and Vgsc is the Galactic escape velocity at the position of 
the Earth). Although the Standard Halo Model, which posits a Maxwell-Boltzmann DM 
velocity distribution, is commonly assumed, it is highly plausible that the DM velocity 
distribution at Earth differs from this, perhaps significantly.' 

This has the consequence that, for elastic scattering, the energy recoil spectrum can 
be anything of the form Fx{q)‘^Fx{q)‘^ multiplied by a non-increasing distribution, as 

® While it is possible that a model with an intermediate-mass mediator would have additional phe¬ 
nomenology distinguishing it from the case of large composite DM, this would not have to be the case. 
Though a light mediator will give rise to self-interactions between galactic DM particles, for heavy (and 
thus dilute) DM these will not be frequent enough to have detectable effects on halo shapes. The direct 
detection event rate is set by the products of the squared couplings for the SM and the dark sector, so by 
making the dark sector coupling large, we could make the required SM coupling very small. Current direct 
detection experiments are sensitive enough that any future signals would imply a minimum value of the 
SM coupling small enough to be significantly below any direct production constraints. Of course, the SM 
coupling may be above this minimum if the dark sector coupling is also small, so direct production signals 
are not ruled out. Also, as discussed in [10], direct bounds on the SM couplings of light (m 100 MeV) 
states generally imply that they cannot lose any initial cosmological energy density they have to the SM 
sufficiently fast (except possibly to neutrinos), requiring the introduction of additional light hidden sector 
states that persist to the present day. 

^Simulations of dark halo formation indicate that non-Maxwell-Boltzmann velocity distributions may 
generically arise [20], and also suggest the possibility of a co-rotating DM disc [21], along with sub-structure 
including cold streams [22]. 
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Figure 3. Left: Example (non-dimensionalised) number distribution of DM states arising from 
nucleosynthesis-type process (from [10]). Right: form factors for DM sizes corresponding to points 
on left-hand plots (upper curve to blue circle, lower curve to red square), and (black) effective form 
factor for whole number distribution. Dashed line shows form factor for light mediator particle, 
F{q) oc l/[q^ + w^), with m = 16MeV. 


discussed in [23]. So, even though the energy recoil spectrum for the dark form factor 
looked different from that for a WIMP assuming the same DM veloeity distribution, by 
changing the velocity distribution we could make a WIMP mimic the (averaged) form 
factor spectrum, since this is decreasing with Er. 

However, if we have data from multiple experiments, each of which uses a different 
type of SM target nucleus, this degeneracy can be lifted. On the assumption that there is 
no dark form factor, the event distributions (once we have compensated for the different 
target nuclei) with respect to Umin should be the same for each of the experiments [23-28] 
— any disagreement indicates the presence of some extra effect. For a dark form factor, 
since g(umin) = ‘^RXNVmin, changing fixN by changing mjg will change the range of Fx{q) 
that we sample (significantly so if the DM mass is larger than those of the SM target 
nuclei). Ref [15] performs a multi-target analysis along these lines for various DM models, 
and a forthcoming paper [29] applies this method to the case of heavy DM coupling via an 
intermediate-mass mediator, which as noted above gives form factors that are very similar 
to those which may be obtained via averaging over Bessel-function form factors. It is 
found® that such a form factor can be distinguished from contact-like scattering with only 
a handful of events in multiple detectors, while effective exposures of only around 1.0 ton 
years in multiple detectors could, assuming DM-SM cross sections close to current upper 
bounds, allow the determination of the mediator mass (which, in the composite case, would 
correspond to the average size of the DM) to around 25% accuracy. 

2.6 Detectability of a rising energy recoil spectrum 

While Sections 2.4 and 2.5 have considered the situation of a distribution over DM sizes, the 
more striking scenario of DM states being concentrated around a single size, giving rise to 
peaks and troughs in the energy recoil spectrum as per Figure 2, has qualitatively different 

®We thank John Cherry for providing plots of the exposure time needed to distinguish our effective form 
factor. 
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Figure 4. Left: Energy recoil spectrum (black) for a Bessel-function form factor corresponding to 
radius 50 fm, with other parameters as per left-hand plot of Figure 2. Blue points (solid lines) are a 
particular sample of 50 events from this distribution, binned with an energy width of 0.3 keV. Red 
points (dashed lines) show the interval of points that is worst fit by a non-increasing distribution (see 
text). The p-value for the test described in Appendix C is 0.005. Right: Cumulative distribution 
functions (CDFs) for p-value from test described in Appendix C, for 30 samples (blue), 50 samples 
(red), and 100 samples (yellow). Grey (dashed) curve shows the CDF for a uniform underlying 
distribution. The upper and lower dashed lines shows the p values corresponding to 2a and 3 ct 
significance for rejecting the hypothesis of a non-increasing distribution. 


detectability in direct detection experiments. As described in the previous section, the 
energy recoil spectrum for elastic scattering can be anything of the form F{q)‘^ multiplied by 
a non-increasing distribution. While this allows non-increasing sections of an energy recoil 
spectrum to be explained by a combination of form factor and DM velocity distribution, 
any rises in the spectrum must come from the form factor (or from the presence of inelastic 
scattering). Though summing over multiple Bessel-function widths will generally smooth 
out the individual peaks into a falling distribution, if the particle sizes were clustered mostly 
around a particular value,® then the troughs/peaks of such an effective form factor could 
give a rising recoil spectrum. Assuming a detector has sufficiently good energy resolution 
so as not to smooth out the troughs in the underlying recoil spectrum, it will be possible, 
with sufficiently many events, to rule out point-like elastic scattering. Furthermore, if 
both the falling and rising parts of a trough were visible, this would be a clear sign of a 
more complicated momentum-dependent form factor, or of some combination of inelastic 
scattering modes. 

We can test whether a given set of events came from a non-increasing distribution 
by locating the energy interval with the ‘worst’ bias towards its high-energy end — here 
measured simply by the average energy of the events — and asking what the probability 
is that a given sample from a candidate non-increasing distribution would have an interval 
that extreme (explained in more detail in Appendix C). Figure 4 displays the results of 
applying this test to simulated data from a model with a Bessel-function dark form factor, 
scaled so that a number of the peaks/troughs are visible in the energy recoil spectrum. As 

® For example, due to the binding energy per constituent reaching a maximum value at some size, and 
decreasing after that (as for iron in the SM), rendering larger states unstable to fission. 
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illustrated, around 12 expected events in the most significant rising section, corresponding 


here to ~ 50 events overall, are sufficient to obtain a two-sigma exclusion of the non¬ 


increasing hypothesis in the majority of cases. While a realistic analysis would need to 
take into account the distribution of background events and other issues, this illustrates 
that relatively few events from only a single detector may suffice to give very interesting 
physical information about DM properties. 

3 Low-energy excitations &; inelastic scattering 
3.1 Properties of low-energy excitations 

Another feature generic to large composite states is the presence of low-lying modes, leading 
to the possibility of inelastic scattering in which these modes are excited. Such modes can 
be broadly characterised as either of “single-particle” or “collective” type. We are most 
interested in the collective modes since, as we will discuss, their excitation amplitudes can 
be coherently enhanced compared to the single-particle excitations, and thus can dominate 
the inelastic scattering rate if they are of low enough energy to be excited. 

There are a number of forms such collective excitations could take. The simplest and 
most generic type of collective excitations are simply vibrational modes of either the bulk 
or of the surface of the large composite state. Alternatively, if the constituents have spins, 
moments, or other intrinsic properties, then collective oscillations of these properties, such 
as spin waves, can occur. In addition, if the ground state is not spherically symmetric, so the 
composite state is deformed in shape as happens for large SM nuclei, then there will be low- 
lying rotational modes as well. For definiteness, in this paper we will focus on the inelastic 
scattering involving the most universal and model-independent of these collective modes, 
namely, the surface and bulk vibrational modes that are of long wavelength compared to 
the scale of the individual constituents. 

For large composite states consisting of roughly uniform ‘nuclear matter’, the lowest- 
energy density waves will generally be surface waves, i.e. volume-preserving oscillations 
of the ‘nuclear surface’ [30]. The reason for this is that the energy of the lowest bulk 
compressional excitations is set by the speed of sound for such waves, Cc ~ where 

K parameterises the compressibility of the nuclear material, and p is the density. Taking 
into account that the lowest possible wavenumber is /c ~ ^/R-, the energy of the low-lying 
compressional excitations is given by 



(3.1) 


On the other hand, the speed of volume-preserving surface capillary waves is set by Cg ~ 
where a is the surface tension (this also holds in relativistic hydrodynamics [31]), so 



the energy of the low-lying surface excitations is given by 



(3.2) 
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Using the facts that R rsj where i?i is the length scale of a single constituent, 

and that both the surface tension and bulk compressibility are set by the same underlying 
interaction strength between the basic constituents so we have a ~ RiK, leads to 


8Ec 


. 


(3.3) 


Thus for large enough A the energy of the bulk compressional modes is well separated from 
the lower-lying surface modes. 

We now summarise the dynamics of the surface modes more quantitatively. Classically, 
for small-amplitude surface waves of a homogeneous, incompressible, sharp-edged fluid 
droplet, we have p{r,6,(j)) constant inside R'{9,(p,t) and zero outside, with 


R'{e,(j),t) = R 


i 

i+EE 

l>2 m=—l 


(3.4) 


for amplitude coefficients The I = 1 modes are removed as they correspond to 

oscillations of the CoM position, while the I = 0 mode is removed as it is just the monopole 
compression oscillation.^*^ The total Hamiltonian of the surface excitations is 

= 2 X/ + Ci\aim{t)f'j , (3.5) 

l>2, 


where Bi = pR?/I and Ci = {I — !)(/ -|- 2)R?‘a are the “mass” and “stiffness” parameters of 
each mode. 

Quantising this system leads to each mode being an independent quantum harmonic 
oscillator. The mode frequencies, of the oscillations are independent of the azimuthal 
spherical harmonic parameter m and are given by [31] 


[Cl _ (l{l + 2){l-l)aV/^ 
y [ R^p ) 


(3.6) 


The overall excitation energy spectrum AE^^i} = ^y the occupation 

numbers n/ of the {21 + l)-fold-degenerate /-modes. It will be important for our discussion 
of the coherent excitation probability of these modes that the simple harmonic oscillator¬ 
like wavefunctions associated to each mode are characterised by a length scale, Rei, where 
the dimensionless parameter ei is given by 


1 

y/2BluJi 


(3.7) 


(For an individual mode with occupation number n; 3> 1 the corresponding classical am¬ 
plitude of oscillation simply scales as \ai\ ~ ei^/nl.) In terms of the underlying parameters 
of the composite state 


ei oc 


^- 7/12 


1 

MiRlP 


1/4 


(3.8) 


**11 the amplitude of the surface waves were not small then volume-preserving oscillations require aoo = 


-E, 


\aim\ , but this constraint can be ignored at leading order. 
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where /3 ~ cj-Rf is the energy associated to surface tension, and Mi is the mass of individual 
constituents. As discussed in [32], the hydrodynamic values of Bi and Ci can differ by 
factors of 0(10) from those obtained experimentally for SM nuclei, due to the effects of shell 
structure etc. However, the hydrodynamic approximation should give a good qualitative 
guide to the properties of the low-lying modes of large composite states. 

Finally, as mentioned above, as well as collective modes there may also be modes 
corresponding to the excitation of ‘single constituents’. If the interior of the state consists of 
degenerate fermionic matter, then excitations which move a state from just below the Fermi 
surface, through a small change in momentum q approximately tangential to the surface, 
will look like quasi-particles of energy q^/(2m*), where the quasi-particle mass m* will often 
be of order the constituent mass [33]. Since the smallest allowed momentum change is g ~ 
1/R, the lowest-energy quasi-particle excitation will have energy ui ~ 2 R^m ^ so 

for large A will be of lower energy than the collective surface modes. However, as discussed 
in the next section, in the regimes of interest in this paper, the cross sections for exciting 
them are much lower. Additionally, the properties of these modes are model-dependent — 
for example, pairing interactions between constituents could lead to superfluidity, as arises 
in forms of SM nuclear matter. 

3.2 Inelastic scattering form factors 

The presence of these low-lying modes means that even low-velocity scattering processes 
may have sufficient energy to be inelastic. In particular, scattering off the long-wavelength 
collective modes, since these involve all of the constituents, will be coherently enhanced as 
for elastic scattering, and have their own momentum-dependent form factors. 

Generally, this form factor can be calculated from the overlap of initial and final 
state wavefunctions. As derived in Appendix A, if the scattering occurs through a scalar 
contact interaction with the constituents of the DM state, then (returning to the notation 
of Section 2.1) the form factor for inelastic scattering into a one-phonon surface mode of 
angular momentum number I is, as given in equation A.8, 

34 

F{q) = ^i^{2l + l)^/hiji{qR) (3.9) 

vdvr 

(as compared to equation 2.2 for elastic scattering), where e/ is the natural amplitude 
associated with oscillations in that mode. This result is correct to first order in e;, and 
valid for qR < 1/e/ (beyond that, the wavefunction overlap can be computed numerically). 
Excitations of multiple phonons are associated with further factors of e/. Comparing elastic 
to inelastic scattering, we see that the latter has form factors corresponding to higher 
spherical Bessel functions, and is suppressed by powers of the natural amplitude of the 
surface modes; as per equation 3.8, this is parametrically small for large A. 

The cross sections for scattering off single-particle excitations add incoherently, giving 
a (X A. There may also be further suppression factors. The case of degenerate fermionic 
matter may be well-approximated by scattering off a Fermi gas of non-interacting quasi¬ 
particles — as discussed in Appendix B, this is, for low momentum transfers, strongly 
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suppressed by degeneracy factors. Generally, incoherent scattering should only be signif¬ 
icant for momentum transfers comparable to the Fermi momentum (and energy transfers 
comparable to the Fermi energy), or for very large composite states where scattering from 
collective modes is highly suppressed (Appendix B). 


3.3 Inelastic recoil spectra 

If collisions are inelastic — say, the DM state is excited from mx to mx + — then the 

minimum velocity [34] for which we can obtain a given recoil energy is 


^min 


1 


/ itinEr 

V ^^XN 



(3.10) 


We have (IEr = where p and p' are the initial and final momentum in the CoM 

frame. The final state phase space available is also proportional to p' , so the differential 
rates for elastic and inelastic scattering are related by 


dai dae |A4ip 
d^ ^ dER |MeP ’ 


(3.11) 


where the A4 are the matrix elements. Aside from differences in the momentum-dependence 
of the matrix elements, giving different form factors, the main qualitative difference from 
elastic scattering is that Umin is no longer monotonically increasing with Er, instead having 
a minimum at Er = dpxN/^N- This is in contrast to elastic scattering, where the fact 
that that Umin is monotonically increasing had important consequences (see Sections 2.5 
and 2.6) for data analysis. 

As per above, we generically expect large composite states to have low-energy exci¬ 
tations, which may be excited by scattering with SM nuclei. The energy recoil spectrum 
will then be a sum over the spectra for scattering into each of these states (including the 
ground state, which gives elastic scattering). Labelling the states by A, 

= X-^f^XnTjv(g)^ ^g(Umin,A(^iR))|v:^-^X,A(9)^, (3-12) 


where ^ gives the ratio of matrix elements (c.f. equation 2.7), and the FA,a have a common 
normalisation such that elastic scattering has F(0) = 1. 

Figure 5 shows an example of the energy recoil spectrum arising from a uniform-density 
composite state of the kind discussed in Section 3.1, in the regime where elastic scattering, 
and inelastic scattering into the first few surface modes, are the dominant effects. Due 
to the relative suppression of inelastic modes, from both form factor effects and the small 
amplitude of surface oscillations, elastic scattering generally dominates in such scenarios. 

Finally we comment that the above analysis assumed that the DM particles are dom¬ 
inantly in the ground state. If de-excitation times are long enough, there is a possibility 

Equation 3.10 is accurate up to multiplicative correction of size at most ~ \J~EJM — this always gives 
a good approximation, since <5 < |/ijv]vUmax is required for any excitations to occur, and galactic escape 
velocity is deeply non-relativistic. 
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Figure 5. Energy recoil spectrum (Germanium target) for spherical DM state with parameters 
as per Figure 2, incorporating elastic scattering and inelastic excitations of surface modes (black). 
Blue (dashed) and red (dot-dashed) curves show contributions from elastic and inelastic scattering 
respectively. The energy of the first surface mode is taken to be 28keV, so a number of modes 
contribute at higher recoil energies, flattening out the inelastic recoil spectrum. The amplitude of 
the inelastic modes is set a factor ~ 10 above the hydrodynamic estimate from Section 3.1, for 
illustrative purposes. Note that the form factor for the first excited mode is 90 degrees out of phase 
from that for elastic scattering — this effect can be seen in inelastic SM-SM scattering, e.g. Figure 
1 of [32]. 



that there may be a significant cosmological population of excited DM states, altering the 
inelastic scattering phenomenology with exothermic interactions being possible. We leave 
discussion of this highly model-dependent possibility for future work. 

4 Astrophysical capture 

In addition to potentially being visible in direct detection experiments, interactions between 
DM and SM matter may also lead to the capture of DM by astrophysical objects, the most 
interesting generally being various kinds of stars. Once enough DM has accumulated inside 
the star, then either self-interactions among the DM (annihilations, self-scattering, etc.), or 
more complicated DM-SM interactions, may in some models lead to observable alterations 
of stellar properties. For definiteness, we consider the case of asymmetric DM, in which 
DM may build up in the star without being destroyed through annihilations. 

For spatially-extended DM states, the possibility of large cross sections for inelastic 
DM-DM interactions may give rise to modifications of the distribution of captured DM 
inside the star, analogously to the possibility of effects on astrophysical halo shapes men¬ 
tioned in [10]. In models with large velocity kicks after exothermic collisions, these could 
take the form of ejecting captured DM from the star, as discussed in e.g. [9]. Alternatively, 
other models could lead to the contraction of the captured DM distribution, potentially 
all of the way down to a very dense configuration at the centre of the star. Since we have 
been considering composite states composed of similar constituents, some kind of repulsive 
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self-interaction countering the attractive binding forces is already required in such models, 
and will generally remove the danger of such a DM configuration collapsing into a black 
hole.^^ In the absence of any further dynamics, there are no obvious externally observable 
consequences of a high-DM-density central region. In particular, the release of binding en¬ 
ergy from fusions will have no significant effect on the star—for symmetric DM, almost all 
of the DM captured by the star can generally annihilate with no observable heating, and we 
are only injecting a small fraction of this energy. However, in the case of run-away fusions, 
there may be inelastic collisions between extremely large composite states, which could 
result in the release of large amounts of energy, in small volumes and on rapid timescales. 
If some of these de-excitation products couple to the SM strongly enough, the SM energy 
injection resulting from these processes may have consequences — for example, [37] points 
out that sufficiently fast and localised energy injection inside a white dwarf may ignite a 
Type la supernova, even in a sub-Chandrasekhar mass dwarf. 


4.1 Effect of dark form factors on capture rate 

The composite nature of the DM may also affect the initial capture rate, through the 
momentum-dependence of scattering cross sections. From [38], the capture rate per unit 
volume is given by 

(4,1) 

dV Jq u 

where u is the velocity at infinity of the DM particle relative to the Sun, f(u) is the DM 
speed distribution at infinity, w is the velocity of the DM particle at the scattering location, 
and i}{w) is the scattering rate to less than escape velocity off the material at that location. 
The scattering rate is 


D('U;) — 'y ^ TliW 0 (Ejiiaxji Emin) 


, ^ dai 


(4.2) 


where n* is the number density of SM nuclear species i, = ]^mxu^ is the minimum 

energy loss required for capture, and Emax.i = with fjLi = niimx/{nii -|- mx) is the 

maximum kinematically-allowed energy loss (taking the ‘cold sun’ approximation where 
the scattering nuclei are at rest). In the case where the scattering involves momentum 
dependent form factors, then using the notation of Section 2.2, 

(T y ■ rQmax.,i 

= / qdqFi{q)^Fx{q)^ ^ (4.3) 


where gmin.i = ^/rmmxu, and g'max.i = 

Combining into a single form factor E{q) = Fi{q)Fx{q), if F((jr) ~ 1 for g < l/R, then 
for qma,y.,iR = ^^iwR < 1 the form factor will have no significant effect. If qmax,iR 1, 
but Q min iR < 1, then assuming that F{q)‘^ falls off faster than the dominant 

may occur in models of very heavy fundamental asymmetric DM [35, 36] 

particular, the Fourier transform of a spherically symmetric distribution with at worst step function 
discontinuities eventually falls off at least as fast as 1/g^. 
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contribution to the integral will be from momenta rsj l/R, giving 


rQma'x,i 

/ qdqF{qf ^ . 


If ^min,z-R ^ 1 as Well, then the integral is 

/'9max,i 


' 9min,i 


qdqFi^q') ^ Qmm,iF(^Qmin,i') -^(i7min,i) iZmaXji 


(4.4) 


(4.5) 


where A will be ~ 1, if qmin,i ^ Qmax,i and F{q)‘^ does not fall off significantly faster 
than a power law in q (if gmin,i > 9 max,i, then A = 0). Turning to the integral over DM 
halo velocities, if the velocity distribution is steeply falling past some typical velocity u and 
approximately constant (in d?u space) below that, as for a Maxwell-Boltzmann distribution, 
then f{u) ~ nxu^/u^i and we can approximate the integral as 


uduwU.{w) . (4.6) 

If Qraax,i{u)R = 2fiiw{u)R < 1 , then no part of this integral has any significant momentum 
dependence (for capture off species i). If qraax,iR A> 1 throughout (i.e. 2niVescR ^ 1, where 
Vesc is the escape velocity), but qmin,i{u)R = ^JnRjnxuR < 1 , then the integral becomes 

uduw^{w)^ni^^'^ . (4.7) 

For momentum-independent scattering with (?max,i 9 min,i(F), 

[ uduwn{w) , (4.8) 

Jo ^ 

(where we assume that Vesc u, so that terms of the form can be ignored), so the 
suppression relative to this case is 4 ^ 2^2 ^^2 = (kR)^’ where k/2 = HiVesc is approximately 
the initial momentum in the CoM frame. 

If Qn ,m ,i(u)R > 1, then 


dC nx r 

Jo 


uduwQ(w) ~ ( / 

2Mf [Jo 


iK^mimxR) I 

udu-—T + 
R^ 


/ udumimxu‘^F{y/minixu)‘^A{y/minixu,2^iw) . (4.9) 

Jl/{yJmimxR) ) 

If F{q)‘^ drops off faster than 1/g^, then this integral is dominated by u l/{y/mmxR)-, 
giving 

axi 1 


(^Xi 


1 


1 


2fsf 2mimxR‘^ 2 yrriimxR'^u^)' ^ 

which has an extra suppression factor. If, as in the Bessel-function case, F{q)‘^ oc 1/q^ for 
large q, then the second integral in equation 4.9 contributes some logarithmic multiple of 
l/{mimxR‘^)- 
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To summarise, for heavy DM {mx rrii) with a radius of R large enough to be 
probed in collisions, kR ^ 1, the elastic capture rate off species i will be suppressed 
relative to momentum-independent scattering by (kR) 2 min (l, , replacing 

the usual kinematic suppression^^ min ^1, ^ (here, C is some logarithmic integral). As 

an example, if we suppose that the Sun (escape velocity 0.002c) has kR ~ 1, then capture 
in white dwarfs {vesc ~ 0.02c) is suppressed by a factor of ~ 10“^, and in neutron stars 
{vesc ~ 0.7c) by ~ 10“®, compared to standard elastic scattering (for white dwarfs, the SM 
nuclear form factor is also important). These suppressions mean that inelastic scattering 
may be the dominant capture process, though we leave the investigation of this to future 
work. 

4.2 Self-interactions of captured dark matter 

One possible effect of large DM self-interaction cross sections is to increase the capture 
rate, in that incoming DM can scatter off already-captured DM as well as off SM nuclei in 
the star. However, for heavy DM, and in locations of reasonable astrophysical DM density, 
it is generally hard to accumulate enough captured DM to have any significant effect on the 
overall capture rate. Supposing that the star has effective capture cross section As for DM 
particles streaming through it, the total number of DM particles it captures in its lifetime 
is Nx ~ AgtgVnx, where nx is the local DM number density and v is the characteristic 
relative velocity (of order galactic orbital velocities, ~ 220 km sec ^). Writing the DM 
self-scattering cross-section as crxx, the total cross sectional area for DM-DM scattering 
is at most 


Ax = (TxxNx ~ {(^xx/mx)tsVpxAs 

^ ( ^xxlmx \ V _ ts _ px 

* \ barn/GeV/ 220 kmsec“^ 5 Gyr 0.3 GeVcm”^ 


(4.11) 

(4.12) 


For elastically-scattering DM with velocity-independent self-scattering cross section, the 
observational limit on the self-scattering cross section is axxlmx ^ Ibarn/GeV (see 
e.g. [39]). That Ax can be of order Ag for this value means that there are possible DM 
models in which self-capture is significant (e.g. [40, 41]). However, for large composite 
DM states, the fact that (Jxxj'^x oc ^x ^ uniform matter means that large states 
generally stand little chance of having significant self-capture. Evaluating (yxxjfnx for a 
composite state of interior density pfe, and constituent number A, 


gaa 0.05barn ^_;^^3 /1 GeV\/I GeVfm 3x^/3 

mA GeV \ mi J \ pb J 


(4.13) 


we see that for SM-like or heavier constituents, large composite states are well below elastic 
self-interaction bounds, so do not give interesting self-capture effects. 

For momentum-independent scattering, if gmin,i(M) > 5max,i, then the integral in equation 4.8 is cut 
off above when — <7max,i, i.e. when u ~ Uesc = Uc (taking pi ~ rrii). This reduces the 

value of the integral by a factor of ~ = rmVgsc/(jnxu^)- 
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However, interactions between already-captured DM states could potentially have in¬ 
teresting effects. Particularly interesting are interactions which may lead to ‘run-away’ 
effects, in which the DM distribution inside the star contracts to a very dense state. This 
could arise either from inelastic collisions which are dissipative (i.e. some of the initial KE 
is lost into de-excitation products), or from fusions, which result in heavier DM states 
that have a correspondingly smaller equilibrium radius inside the star. In Sections 4.2.1 
and 4.2.2, we will perform some approximate calculations to demonstrate the feasibility of 
these scenarios. 


4.2.1 Dissipative collisions 

For dissipative collisions, the dynamics are governed by the rate at which DM-DM collisions 
dissipate energy, versus the rate at which DM-SM collisions add it, re-thermalising the DM. 
As an estimate, if some fraction a of integrated phase space density for DM lies within a 
spatial radius r, then the rate of self-interactions is set by 


Fxx ~ crxxnxvx ~ 


crxx 

J.2 


aNx 


3 

dvr 



r\j 


10 day ^ 



(4.14) 


where we have taken v'j^ ~ , and Ax is as per equation 4.12. In thermal equilibrium 

with the SM matter in the star, the DM states would have an isothermal distribution, 


Px 






r* = 


('47rGmxp* 


1/2 


(4.15) 


where T* and /j* are the temperature and density of the stellar core. Evaluating Txx for 
r = r^:, and putting in values appropriate to the Sun, 


Fxx 


6 yr 


( CFXnImx \ 
VlO-8 pb/TeVy 


mx \2/3 / I GeV fm ^ 

TWJ V Jb 



(4.16) 


where we have used equation D.l for the solar capture rate, and we have replaced axx 
by the geometrical cross section between two spheres of mass mx and internal density 
This interaction timescale is much less than the lifetime of the Sun. So, as long 
as scattering with SM particles does not counteract the increase in the local DM number 
density that comes about from losing energy in inelastic scatterings, it is inconsistent for 
a large fraction of the captured DM to have a basically isotropic steady state phase space 
distribution with most of the density within a sphere of isothermal or smaller radius. Thus, 
if the DM-SM cross section is large enough, the DM distribution may contract down to a 
very dense configuration at the centre of the star. 

It remains to check that SM scatterings do not re-thermalise DM fast enough to 
avoid this contraction. The viscous drag force on a DM state is approximately fdrag ~ 
—mxnxvxcrxN'vx = —I'fnx'^x-, where mx, nx and vx are the mass, number density and 

1 TeV mass sphere of this density has radius ~ 6fm, which is of the order of the inverse momentum 
transfer in the initial DM-nucleus collisions (for solar capture), so there will not be significant momentum 
suppression of the capture rate. 
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average speed of the SM scatterers, axN is the momentum transfer cross section between X 
and the scatterers, and vx is the velocity of X in the thermal rest frame of the scatterers. 
Thus, 7 gives the damping rate at which scatterings return vx to the thermal distribution. 
For parameters relevant to the Solar core. 


mxnxVNcrxN „ i„_4 _i 

7 ~-~ 3 X 10 yr 

mx 


( CFXnlmx \ 
VlO-8 pb/TeVy 


(4.17) 


Thus, as long as each inelastic collision does not dissipate too small a fraction of the initial 
KE (~ lO”'^ for the parameters in equations 4.16 and 4.17), SM-DM thermalisation will be 
too slow to prevent contraction. Also, since the thermalisation timescale is much less than 
the lifetime of the star, we do expect to contract down to isothermal densities in the hrst 
place (in particular, this means that we will lose almost all of any net angular momentum 
the captured DM distribution might have had initially, meaning that rotational support 
of the collapsing distribution will not be a worry). The same separation between self¬ 
interaction rates and thermalisation times may, for suitable parameter ranges, apply in the 
case of white dwarfs and neutron stars as well (however, for neutron stars there is an issue 
of whether thermalisation timescales are less than the lifetime of the star [42]). 


4.2.2 Fusions 

An alternative way to realise run-away contraction is for the DM states to progressively 
increase in mass through fusions, and correspondingly for their isothermal radius within the 
star to decrease. In contrast to the previous section, this scenario relies on thermalisation 
through SM-DM scattering being fast enough. It also relies on fusion cross sections being 
large enough, and on fusions remaining energetically favourable up to very large sizes — 
these features may naturally be realised in models of composite DM such as those discussed 
in [10]. 

As discussed in Appendix B, we expect the coherent momentum transfer cross section 
between DM and SM states to scale as a an oc for states A with pRa ^ 1, where p is 
the characteristic momentum of the SM scatterers. Thus, 'ja — (Ap/Ay^^jAp for A » Ap, 
where Ap is the size for which the DM states are of radius 1/p, i.e. pRap 1, with 7 
given by equation 4.17. In terms of the size Aj of states which have radius ~ 1/pi, where 
Pi is the characteristic momentum of SM scatterers in the initial collision (which occurs 
at around the escape velocity), thermalisation is naively effective (assuming small enough 

injections of kinetic energy on the DM side) for A < 10^®Aj ^ nj-s pb^xeV ) = The 

isothermal distribution for such heavy states may already have volume comparable to their 
saturated volume. If not, then further fusions will generally be fast enough to combine 
most of the DM into a few very large composites — the initial number density, and so rate 
of fusions, will be very high, with build-up continuing as Taa oc A“^A^/^ ~ A“^/^ from 
number density and cross section factors respectively. 

For white dwarfs and neutron stars, we have a similar conclusion that fusions are 
potentially fast enough to form a very dense configuration. In these cases, the star is 
compact enough that even a DM distribution of radius comparable to the entire star can 
still have a high enough self-interaction rate for fusions to combine the majority of the 
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DM into a dense state (though we would still need to worry about dissipating energy and 
angular momentum, especially in the neutron star case). 

We note that the estimates in Section 4.2.1 and 4.2.2 should be viewed as rough plau¬ 
sibility estimates — proper investigation of these issues would require realistic modelling 
of the DM-DM collisions, and of the phase space distribution of DM at each stage of the 
process. In particular, this section has ignored the possibility that DM collision types other 
than fusions are important (e.g. fragmentations etc.), and also assumed that the veloc¬ 
ity kicks imparted by fusion de-excitation products are small enough to be re-thermalised 
quickly. However, the point was merely to illustrate that run-away contraction is a plausible 
possibility. 

5 Conclusions 

In this paper, we have investigated some of the consequences that follow if a proportion of 
DM is composed of large composite states — specifically, states consisting of a large number 
of constituents forming an extended, semi-uniform object with a saturated density. 

The spatial extension of these objects introduces a dark form factor into elastic scatter¬ 
ing amplitudes as discussed in Section 2. This has two effects: First, coherently enhanced 
scattering rates with the result that collider bounds on DM cross sections are effectively 
weakened by factors of up to 1/H (Section 2.3). Second, the introduction of an additional 
characteristic momentum dependence, whose effects on direct detection energy recoil spec¬ 
tra, for DM states with radii large compared to SM nuclei, are considered in Sections 2.3 
and 2.4. We find that, in the most visible cases, such signals may be distinguishable from 
elastic scattering, in a halo-velocity-independent manner, with only 0(50) scattering events 
(Section 2.6). 

Large composite states will generically have long-wavelength collective excitations, 
which give rise to the possibility of low-energy, coherently enhanced inelastic scattering 
processes, as discussed in Section 3. Though such processes will generally be sub-dominant 
in direct detection experiments (Section 3.3), in specific models they may have signatures, 
in direct detection or astrophysically. 

We considered the capture of DM by stars in Section 4. The presence of dark form 
factors and inelastic scattering may alter the capture rate. In addition we investigated 
the possibility of inelastic interactions between captured DM states having an effect on 
the DM distribution inside the star. Strikingly, for plausible types of self-interactions, we 
argue that it is possible for almost all of the captured DM to accumulate into a very dense 
conhguration (Sections 4.2.1 and 4.2.2). 

Finally, we emphasise that our calculations have focused on what we believe to be the 
leading, most model-independent features associated with the interaction of large composite 
DM states with the SM sector. In specific realisations of large composite DM there may 
well be other even more striking signals. It would clearly be interesting to construct specific 
models which realise large composite DM and investigate the resulting phenomenology. 
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A Matrix elements for inelastic scattering 


Suppose that we have a scattering process in which a state i, with initial momentum k, 
scatters off a state I with initial momentum p (all non-relativistic), resulting in final states 
/ and F with momenta k' and p' respectively. The matrix element for this, in the Born 
approximation, is 


A = (/,k';F,p'|iLint|i,k;/,p) 

= (/,k;T,p|A‘i-("^-"^)iLi„t|i,k;/,p) 
= [ drF^■’^Ves{r,^r) , 


(A.l) 

(A.2) 

(A.3) 


where k' — k = p — p^ = q is the momentum transfer, v is the relative velocity of i and I, 
and (xj; xp’|FIint|k; p) = (27r)“^I4fr(x/ — x^?, That is, the matrix element is 

given, as a function of the momentum transfer q, by the Fourier transform of a (possibly 
velocity-dependent) effective potential Ves, whose form will depend on the interaction 
Hamiltonian. If the initial or final states have directional properties (polarisations), I4fr 
may also depend on those. 

As discussed in Section 3.1, the surface modes of an incompressible liquid drop can be 
obtained by quantising the classical surface oscillations, in which the surface is displaced 
as in equation 3.4. Classically, for a scalar interaction between a plane-wave scatterer and 
the nuclear matter, the Fourier transform of the interaction potential is given by 


F{aim)= ! d^vF^-^Vo 




= Vo dn dr '^{21' + ji^{qr)Pi:(cos 9) , 


(A.4) 

(A.5) 


where we have taken q to be in the z direction. For m ^ 0, this integral is clearly zero. 
Expanding to first order in a given aio, and eliding Vq, 

F{aio) = F{0) + V^ [ dQaioR^Yio{9) ^(2/' + ji,{qR)YMe) (A.6) 


3A 


= F{0)+aio^{2l + l)Fhljl{qR) . 
vdvr 


(A.7) 
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Treating each mode as a harmonic oscillator, we have d/o = ^i{^i + where ei is the 
dimensionless amplitude from equation 3.7, and a| is the creation operator for the 1,0 
mode. So the matrix element between the ground state and the first excited state is, to 
first order in e;, given by 

3/1 

(lHF(d)|0) = ^{21 + . (A.8) 

Vd-TT 

Higher phonon number states are associated with further factors of ei (for the appropriate 
I numbers), and more factors of qR in front of the spherical Bessel function. For the 
expansion to make sense, we must have that ji{qr) varies slowly over the interval ii(l ± e), 
so we need qR < 7r/e;. Beyond this approximation, 

(lz|F(d)|0) = J da'ijjl{a)'tlJo{a)F{a) , (A.9) 

where the V’n are the oscillator wavefunctions (for a harmonic oscillator, Hermite polynomi¬ 
als multiplied by an exponential), and F{a) is from equation A.5 (note that the oscillator 
properties derived in Section 3.1 are also to first order in ei). In particular, comparing 
equation A.8 to the elastic scattering form factor 2>Aji{qR)/{qR) from equation 2.2, both 
are ~ Aef at qR ~ 1/e/, but the correct form of the inelastic form factor from equation A.9 
drops off faster than the elastic form factor beyond this, since the wavefunction for the 
surface displacement integrates over multiple, cancelling, Bessel-function periods. 


B High-momentum scattering from large composite states 

In this Appendix, we consider the scattering behaviour of large composite states when the 
momentum p of the incoming scatterer is such that pR ^ 1, where R is the radius of the 
composite state. 


B.l Coherent scattering 

The elastic scattering form factor should fall off at least as fast as {qR)~^. For thermalisation- 
type processes, we are most interested in the rate at which SM scatterings exchange energy- 
momentum with the composite state, which for dominantly soft scattering will be set 
parametrically by the momentum transfer cross section an = f (1 — cos0)^dQ. Writing 
m ~ using q^ = 2p?v‘^{l — cosO) in the CoM frame. 


1 

O'tr — 2*^0 


d COS 0(1 


cos 9) 


fjqR) 

{qRY 


1 ctq 

2 {pvR)'^ 



dq 

q 


f{qR) ■ 


The dimensionless integral on the RHS provides a logarithmic factor A, giving 


(B.l) 


(B.2) 


Since (Tq oc A^, we have an ~ for A large enough that kR ^ 1 (with A changing only 
logarithmically with kR), as used in Section 4.2.2. 
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The above assumed that the form factor was that appropriate for plane-wave scattering. 
If both the composite state and the SM scatterer are better-localised than the size of the 
composite state, then the wavefunction overlap in the form factor will not probe the full 
composite state, giving different results. 

For inelastic scattering, taking the composite state to be much heavier and much larger 
than the scattering state, we can approximate the scattering as being against an infinite 
uniform medium (assuming that p is small enough that it does not resolve structure on 
the scale of individual constituents). Since the collective modes are linearly dispersing, 
Ld = Cck, then by the usual pseudo-momentum and energy conservation considerations, 
scatterers can only excite these if their velocity relative to the medium is greater than Cc 
(as per superfluids). Putting in some illustrative numbers, the speed of CNO nuclei inside 
the Sun is u ~ 3 x lO^'^c, while the speed of sound for compressional modes of SM-like 
nuclear matter is ~ 0.1c. 

For a composite state of finite radius, this corresponds to the fact that the modes of 
small enough energy to be excited have k q, so the wavefunction overlap in the form 
factor is very small (the q wave oscillates much faster than the k wave, giving a large 
cancellation). With collective surface modes, as explained in Appendix A, for eiqR S> 1 
there is again a large cancellation. 

The above comments assume that the composite state is in its ground state. If there 
is some non-zero occupation number for high-wavenumber modes, then down-scatterings 
of these are energetically permitted, and will not suffer from the cancellation suppression 
described. We do not consider the case of scattering against excited states here, assum¬ 
ing that the de-excitation times are much shorter than the times between composite-SM 
collisions (for this to be the case, de-excitation will generally have to be to hidden sector 
states). 

B.2 Incoherent scattering 

As discussed in Section 3.1, as well as the collective modes, there may also be low-lying 
modes corresponding to the excitation of ‘single constituents’ — in the case of degenerate 
fermionic matter, excitations in which particles just below the Fermi surface are scattered 
to just above it. By the usual Fermi-liquid theory, ignoring interactions and approximating 
the scattering as occurring from a non-interacting Fermi gas should provide a good first 
approximation. [42] discusses scattering from a degenerate Fermi gas, finding that the 
scattering rate for low-energy scatterers is independent of the Fermi momentum (as can 
be seen from geometrical considerations). This can result in large suppression factors 
compared to the naive scattering rate given by anv — for scatterers with momentum k <C 
Pf, and also low velocity compared to the Fermi velocity, we have an effective suppression 
(for energy transfer rates) of ~ (, where mx is the mass of the scatterer, mn is 
the effective mass of the quasi-particle, and p = mxfnn/{mx + mn) is the reduced mass. 

Since the cross sections for scattering off single-particle excitations also add incoher¬ 
ently, a oc A, coherent scattering should dominate unless the size of the state is very large, 
so that as reviewed in the previous sub-section, the coherent transfer cross section grows 
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slower than A. Making an estimate for SM-like nuclear matter (individual constituents of 
mass mi = 1 GeV, bulk nuclear matter density p = 1 GeV/ fm^) scattering off solar mate¬ 
rial, we would need, very roughly, A > 10^^ for incoherent scattering to start dominating 
the transfer cross section (corresponding to i? > 40 pm). This is well above the size be¬ 
yond which thermalisations were found to be ineffective in Section 4.2.2. However, it does 
illustrate why, for example, scattering of DM in neutron stars is dominated by incoherent 
scattering. 

C Statistical identification of rising distributions 

Suppose that we have some (one-dimensional) data points, which we assume are IID sam¬ 
ples from some probability distribution. We wish to test the hypothesis that this distribu¬ 
tion is non-increasing, with respect to some appropriate function of the parameter. 

If we expect plausible alternative distributions to feature only one prominent rising 
segment (for example, the recoil spectra considered in Section 2 have successive Bessel 
function peaks suppressed by both the natural fall-off of the Bessel function, and the 
velocity distribution), a sensible approach is to locate the interval with the ‘worst’ bias 
towards its right-hand end, and ask what the probability is that a given sample from a 
candidate non-increasing distribution would have an interval ‘that extreme’.^® 

A simple measure of how biased an interval is towards its right-hand end is simply 
the average position of the points within it. The non-increasing distribution maximising 
the probability of right-biased points is clearly the uniform distribution on the interval. If 
observations are binned, we can find the exact distribution for the average of bin mid-point 
positions, for some number n of samples, by performing the convolution of the (binned) 
uniform distribution on the interval with itself n times. In particular, if the bins are of 
uniform widths, this distribution can be computed analytically. If observations are not 
binned, then the distribution for the average of the positions is the (rescaled) Irwin-Hall 
distribution for n points. In both the discrete and continuous cases, the null distribution 
is approximately normal for large n, with variance l/\/l2n (on the interval [0,1]). 

So, if we have n samples with a mean position of x (rescaling the interval to be of 
width 1), we can use the appropriate distribution to find the probability of a mean position 
> X arising from a uniform distribution on that interval. The test statistic for the whole 
sample is then the minimum such p-value for each sub-interval within the sample (since it 
is clear that the worst sub-interval will always terminate either at the left sample end, or at 
the location of a point, there are a finite number of sub-intervals to test). To determine the 
distribution of this test statistic under a candidate non-increasing distribution, Monte Carlo 
simulation can be used. A uniform distribution is clearly the non-increasing distribution 
most likely to produce fake rises, and while the data may indicate that the underlying 
distribution is far from uniform, for reasonably small numbers of samples the p-values 
obtained by adopting a uniform null distribution are only slightly worse than those obtained 

^®This would be a poor approach if we expected the rises in the distribution to be e.g. a small periodic 
signal super-imposed on some larger background, or more generally any small but structured deviation from 
a larger background. 
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by allowing a free null distribution constrained by the fit to data. In particular, we adopt 
a uniform null distribution for the calculations in Figure 4, and find that the p-values are 
generally no more than a factor ~ 2 worse. 

One issue is that, in many physical cases, there will be some resolution associated 
with our points. Then, the distribution our samples are drawn from is some (positive) 
underlying distribution convolved by a (positive) detector response function, so must be 
smooth on scales of order the resolution. This means that rightwards bunching on such 
scales must be spurious, so should not be considered in our test. The simplest way to solve 
this issue, if we expect rising features in plausible alternative distributions to be on scales 
larger than the resolution, is to bin the points on around the resolution scale. 

Quantitatively, the variance for the average of n uniform random variables on [0,1] is 
and the expected average position from a linearly rising distribution is | + g, so we 
would expect to exclude the uniform distribution at a significance of ~ 0.6y^ sigma. In 
the example from Figure 4, with 100 events overall, we expect around 26 in the first rise, 
so we would expect to obtain around a 3a exclusion, assuming that the worst intervals 
arising from the null distribution are of approximately that many events. This is indeed 
what we observe, with the p-value CDF in the 100-sample case following very closely this 
normal approximation (the yellow curve in Figure 4, lying almost entirely underneath the 
yellow points). The CDFs for the 30 and 50 sample cases are not so well approximated, 
since the ‘worst interval’ is more variable with small numbers of events. 

D Solar capture of heavy WIMPS 

Here, we give a very brief review of the capture of weakly-interacting, massive dark matter 
particles by the Sun. For X particles much heavier than any of the relevant SM nuclei 
within the Sun (here, this means mx ^ 400 GeV), the rate of capture by the Sun scales as 
Cq ~ axnPx'm-x^, where axn is the scattering cross section with SM nucleons (assumed to 
be elastic and spin-independent) [38, 43].^^ Summing over capture rates from the various 
elements in the Sun (the main contributions are from the CNO elements) using a standard 
solar model [44], we obtain 

Cs = 3xl0«sec-‘(j^)(i^) , (D.l) 

(taking the couplings to protons and neutrons to be the same). We have assumed that the 
DM velocity distribution follows the Standard Halo Model—since, for high-mass WIMPs, 
only the low-velocity part of the distribution can be captured by scattering events, modifi¬ 
cations that affect the low-velocity distribution will alter the capture rate, direct detection 
experiments imply that axn ^ 10~^ [45], and the age of the Sun is t© ~ 5Gyr, so 

One factor of comes from the number density nx = pxm^^, while another arises from the fact 
that heavier particles lose a smaller fraction of their kinetic energy in collisions with SM nuclei, so only 
the low-speed part of the WIMP velocity distribution can be captured. This also assumes that axn is low 
enough that the Sun is optically thin to X particles, which is the case for axn < 10“'^ pb. 
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the number of X particles captured is < 5 x 10^® (^m^) ■ ^ fraction of the total flux 

of X particles hitting the Sun, we capture ~ 3 x 10“^ ( lo-^pb^Tev ) • 

Here, ‘captured’ means that the X particles are in gravitationally bound orbits passing 
through the Sun. Subsequent scatterings with material in the Sun will reduce the size of 
these orbits further, and eventually the X particles will (ignoring other interactions) settle 
into an isothermal distribution px ~ with [46] 


r* = 


/-37©-I 

\A'KGmxPQ J 


~ 2 X lO“^i?0 


mx 




107 K 


^150gem 


PQ 


(D.2) 


(where the temperature and density are appropriate to the solar core). If mx 100 TeV, 
then the initial orbits will be small enough that planetary perturbations can be mostly 
neglected, and then thermalisation occurs in less than the lifetime of the Sun if axn ^ 
3 X 10“^^ pb {mxI TeV)^^^ [47]. For larger mx, most of the initial orbits are large enough 
that the effect of Jupiter perturbs them so that they no longer pass through the Sun, 
resulting in most of them never thermalising [47]. 

Finally we remark that the above calculations apply to the case of capture by elastic 
scattering with SM nuclei. If there are sufficiently low-lying excitations that inelastic 
scatterings are possible, in either or both direct detection experiments and solar capture, 
this may change the possible parameter space. (See [48, 49] for related investigations of 
solar capture of inelastic DM.) 
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